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Abstract 

The grand canonical ensemble formulation of the van der Waals equation of state that includes the 
effects of Bose statistics is applied to the equilibrium system of interacting pions. If the attractive 
interaction between pions is large enough, a limiting temperature Tq emerges, i.e., no thermodynamical 
equilibrium is possible at T > Tq. The system pressure p, particle number density n, and energy density 
e remain finite at T = Tq, whereas for T near Tq both the specific heat C = de/dT and the scaled 
variance of particle number fluctuations a;[A] are proportional to {Tq — T)~^D and, thus, go to infinity 
at T Tq. The limiting temperature also corresponds to the softest point of the equation of state, 
i.e., the speed of sound squared cf = dp/de goes to zero as (Tq — Very similar thermodynamical 

behavior takes place in the Hagedorn model for the special choice of a power, namely m~^, in the 
pre-exponential factor of the mass spectrum p{m). 

PACS numbers: 12.40.Ee, 12.40.-y 
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I. INTRODUCTION 


The van der Waals (VDW) equation of state (EoS) is a simple analytical model of the 
pressure function p for the system of particles with both attractive and repulsive interactions. 
In the canonical ensemble (CE), where independent variables are temperature T, volume V, 
and number of particles N, the VDW EoS has the most simple and transparent form, (see, e.g.. 
Refs. mm, 


p(T, n) 


NT 

V-bN ~ 


riT 
1 — hn 


an^ , 


( 1 ) 


where a > 0 and 6 > 0 are the VDW parameters that describe attractive and repulsive in¬ 
teractions, respectively, and n = N/V is the particle number density. The first term in the 
right-hand-side of Eq. ([^ contains the excluded volume correction (for instance, b = 167rr^/3 
with r being the particle hard-core radius), the second term comes from the mean field descrip¬ 
tion of the attractive interactions. 

In order to apply the VDW EoS to systems with variable number of particles the grand 
canonical ensemble (GCE) formulation is needed. This was done for the VDW equation ([^ in 
our recent paper [3]. The GCE formulation of EoS in the form p = p{T, n), including the VDW 
equation, can also be conveniently treated within the thermodynamic mean-field approach 
E]. As the next step, we proposed in Ref. the generalization of the VDW EoS that includes 
effects of the quantum statistics. The nuclear matter was considered in [7j as the system of 
interacting nucleons with Fermi statistics and the VDW EoS. With this model the curve of the 
first order phase transition and its end point (the so called critical point) were obtained. In 
Ref. [8] the fluctuations in the vicinity of the critical point were investigated. 

In the present paper we apply the VDW EoS with quantum statistics to a description of 
the gas of interacting pions. In this case the GCE formulation is really necessary: the number 
of pions is not a conserved quantity and cannot be considered as an independent variable of 
the CE. Note, the EoS of interacting pions was studied during last decades within different 
theoretical approaches (see, e.g., PR3] and references therein). 

The paper is organized as follows. In Sec. [IT] we present the VDW EoS that includes the 


effects of quantum statistics. In Sec. Ill this formalism is applied to the pion gas. The effects 
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related to the limiting temperature are investigated in Sec. IV and a comparison with the 


Hagedorn model is conducted in Sec. M A summary in Sec. VI closes the article. 


II. VDW EQUATION OF STATE FOR QUANTUM STATISTICS 


The pressure and the particle number density for the quantum VDW EoS in the GCE are 
dehned as the following [7j: 

p(r,/i) = - an^{T,^i) , n{T, p) = 

where 

p* = p — bp{T, fi) — abn'^{T, fi) + 2 an{T, fi) , (3) 


with p being the chemical potential which regulates the particle number density in the GCE. 
Quantities and correspond to the quantum ideal gas pressure and particle density, 
respectively: 


p^%T,p*) = 


(Pk 


(27r)3 


exp 


y/m'^ + IP — p* 


-1 


T 


+ p 


= g 


(Pk 

(27r)3 


exp 


^/np + \P — p* 


-1 


V 


( 4 ) 

( 5 ) 


where g is the degeneracy factor, m is the particle mass, p = +1 for Fermi statistics, p = —1 
for Bose statistics. 

The quantum formulation (|^-(|^ of the VDW EoS in the GCE has the following basic 
properties: (i) it transforms to the ideal quantum gas if both a = 0 and 6 = 0; (ii) it becomes 
equivalent to the classical VDW EoS Q in a region of thermodynamical parameters where 
quantum statistics can be neglected; (hi) the entropy found from Eqs. ([^-([^ is a non-negative 
quantity and it goes to zero at T —)■ 0 in accordance with the third law of thermodynamics (see 
Ref. [7j for details). Note, in the Boltzmann approximation, i.e., p = 0 in Eqs. (|^ and ([^, 
EOS obtained in the GCE from (|^-(|^ as p = p(T, n) becomes identical to the VDW equation 
of state given in Eq. ([^. 
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III. PION GAS WITH THE VDW EOS 


In a case of the pion gas one should use the Bose statistics, rj = —Im. Eqs.Q and (|^. The 
pion system with zero value of total electric charge is considered, thus, tt’*', tt”, and 7r° all have 
zero chemical potentials, i.e., /i = 0 in Eq. (|^. As a result, p and n in Eq. (|^ become the 
functions of one variable, T, and Eqs. ([^-(|^ take the form: 


p* 


hp{T) — abn‘^{T) + 2an{T) , 


1 

3 


k^dk 


k^dk 


k^ 


\/ml + 


exp 


ij* 


exp 


+ — p* 

T 


T 


- 1 


-1 



( 6 ) 

(7) 

( 8 ) 


where we set g = 3 for the pion degeneracy factor and = 138 MeV for the pion mass. In 
what follows we take a hxed value of 6 = 0.45 fm^ that corresponds to r = 0.3 fm“^ for the 
pion hard-core radius [HIIS] The value of a is considered as a free model parameter. Note 
also that the entropy density s{T) and energy density s{T) can be calculated from the function 
p(T) as the following: 


s = 


dp 

llT ’ 


e{T) = T 


dp 

df 


- P ■ 


(9) 


We consider T < 160 MeV as a region, where the pion gas may exist. At higher temperatures 
the effects related to the deconhnement are expected to play a major role. For every T, at 
T < 160 MeV, we hnd p* = p*{T) by solving numerically the transcendental equation (|^, the 
functions n{T) and p{T) are then found from Eq. ([^. If Eq. ([^ has more then one solution for 
a given T, only the solution with a larger value of p{T) should be considered as a stable one 
according to the Gibbs criteria. For the VDW EoS in the GCE these criteria were considered 
in Ref. [7] 

The results for p*, p, n, and e as functions of T are shown in Figs. (a)-(d) by solid lines 
for different values of the parameter a. Dotted lines in these hgures correspond to the ideal 


^ In the next section we discuss how the model results depend on the numerical value of b, particularly at 
b^O. 
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Figure 1: The thermodynamical quantities of the pion gas with the VDW EoS as functions of temper¬ 
ature: shifted chemical potential n* (a), particle number density (b), pressure (c), and energy density 
(d). The excluded volume parameter is fixed as 6 = 0.45 fm^. The solid lines correspond to the 
solutions with a/b = 0, 200, 400, 500, 600 MeV (from bottom to top). The dotted lines show the 
behavior of the ideal pion gas, i.e., for both a = 0 and 5 = 0. The points at the limiting temperature 
are depicted by the cross symbols x. In the panel (a) the unstable solutions are depicted by dashed 
line segments with endpoints at T = T' depicted by the filled circles. 

pion gas (i.e., both a = 0 and 5 = 0). In the region of T < 100 MeV (not shown in Fig. 
the interaction effects for the considered a and 5 values are small, i.e., all lines are very close 
to those of the ideal gas. This is due to the fact that the VDW interactions depend on the 
particle number density, and they become negligible at n —)■ 0. 

The repulsive interactions suppress the VDW thermodynamic functions n{T), p{T), and 
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e{T) in comparison to those of the ideal pion gas. This is clearly seen for lines with a = 0 in 
Figs. (b)-(d). The VDW repulsion also introduces the upper limit for the particle number 
density, namely, n < 1/b = 2.21 fm“^. On the other hand, the presence of the attractive 
interactions enhances n{T), p{T), and s{T) in comparison to the ideal pion gas. One observes 
that the VDW EoS with b = 0.45 fm^ and a/b = 200 MeV appears to be close to that of the 
ideal gas, i.e., the actions of repulsive and attractive interactions almost cancel each other out 
in all thermodynamical functions. 

We note that the VDW excluded volume correction is consistent with the virial expansion 
for hard spheres only at low enough densities. At large values of the so-called packing fraction 
Tj = (bn)/A (a fraction of the total volume occupied by the particles of finite size), namely 
7] > 0.1, the VDW equation deviates signihcantly from the equation of state for hard spheres. 
The highest value of the pion number density obtained in the present work is Umax — 0.30 fm“^, 
as seen in Fig. [Tj Therefore, the maximum value of the packing fraction is around r^max = 0.034. 
At such values of p the VDW equation is still fully consistent with the virial expansion for the 
hard spheres (see, e.g.. Fig. 1 in Ref. m) and thus can be used. 


IV. LIMITING TEMPERATURE 

At large values of the parameter a a new phenomenon takes place: a limiting temperature Tq 
emerges, i.e., Eq. (j^ has no stable solutions at T > Tq At the same time in the temperature 
interval T' <T < Tq the unstable solutions for the VDW thermodynamical functions appear. 
Here the temperature T' is determined as a starting point of the Bose-Einstein condensation, 
i.e., p*{T') = rriT^. The values of the thermodynamical functions at the limiting temperature Tq 
are depicted by the cross symbols x in Figs, [^(a)-(d), while the values at T = T' are depicted 
in Fig. (a) by circles, and the unstable solutions are shown by the dashed segments of lines. 
The stable solutions shown in all hgures by solid lines correspond to the Gibbs criteria, i.e., they 
have pressure larger than that of the unstable solutions at equal temperatures. The unstable 
solutions are shown only for a shifted chemical potential /i*(T), in the panel (a) of Fig. [^ 

^ Note that only temperatures T < 160 MeV are considered. 
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Figure 2: (a) The heat capacity, (b) The speed of sound squared (the dotted horizontal line corre¬ 

sponds to = 1/3, obtained for the ideal gas of massless particles), (c) The scaled variance uj[N] 


(10). The solid lines correspond to solutions with b = 0.45 fm^ and a/b = 0, 200, 400, 500, 600 MeV 


from bottom to top in (a) and (c), and from top to bottom in (b). (d) A dependence of the limiting 
temperature Tq on the parameter a for 6 = 0 (dashed line), 6 = 0.45 fm^ (solid line), and 6 = 1 fm^ 
(dashed-dotted line). 


To elucidate a physical origin of the limiting temperature Tq for the pion gas with the VDW 
EoS we calculate the specihc heat, C = de/dT, the speed of sound squared, = dp/de, and 
the scaled variance of particle number fluctuations. 




2an(T, p*) 


(1 — bnY 


a;id(T, p* 


-1 


( 10 ) 
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where 


uJid{T,iJ*) = 


n'^{T,ij,*) 2 tt‘^ Jq 


k^dk 


exp 


+ k‘^ — jj,* 

T 


-2 


( 11 ) 


is the scaled variance of the particle number fluctuations in the ideal pion gas with the chemical 
potential fi*. Equations (10) and ( [II| were obtained for the quantum VDW EoS in Ref. |8]. 
As seen from Figs. (a)-(c), —)■ 0, while both C —)■ oo and Ci;[A^] —)■ cxd at T —)■ Tq. The 

straightforward calculations reveal the power-law behavior at T Tq: both C and a;[A^] are 
proportional to (Tq — and oc (Tq — T)^A_ Fluctuations of the total energy E behave 

very similar to fluctuations of the number of pious N. The limiting temperature corresponds, 
therefore, to the softest point of the EoS (the speed of sound equals zero) and to inhnitely large 
values of the heat capacity and fluctuations. In Fig. (d) a dependence of Tq on parameter a 
is shown for different values of repulsive parameter b: b = 0, b = 0.45 fm^ and 6 = 1 fm^. It is 
seen from Fig. (d) that, at each value of VDW repulsion parameter 6, there is a minimal value 
of VDW attraction parameter a = at such that the limiting temperature T = Tq < 160 MeV 
appears at all a > a^, and Tq decreases with increasing a. At 6 = 0.45 fm^, one hnds that 
To = 140 MeV and 155 MeV at a = 270 MeV fm^ and 225 MeV fm^, respectively. For the 
smaller values of 6, the same values of the limiting temperature, Tq = 140 MeV and 155 MeV, 
are reached at smaller values of a. Particularly, for 6 = 0, the values of Tq = 140 MeV and 
155 MeV are obtained at a = 219 MeV fm^ and 168 MeV fm^, respectively. Note that additional 
information is needed to have further restrictions on the parameters a and 6. From the present 
analysis the most interesting physical behavior with Tq = 150 ± 10 MeV takes place in the 
range of 6 from 0 to 0.45 fm^ and of a from 150 MeV fm^ to 300 MeV fm^. 

The VDW parameters a and 6 were obtained in Ref. [7j for a system of nucleons. Based 
on the ground state properties of the nuclear matter it has been estimated that for nucleons 
6 ~ 3.42 fm^ and a ~ 329 MeV fm^. Therefore, the numerical values of both 6 and a VDW 
parameters in the pion gas considered in the present study are smaller than those in the system 
of interacting nucleons. 

In principle, the information about the VDW parameters a and 6 could be inferred from 
the TTTr-scattering data. Note that second cluster integral for the equilibrium pion gas was 
calculated in Ref. [IB]. The phase shifts in ttyt scattering were described by assuming that both 








the hard-core repulsion and resonance attraction are equally important and should be taken 


into account simultaneously. From the specihc calculation in Ref. [16] it could be inferred that 


parameter a may take values up to few hundred MeV fm^, and could in general be temperature 
dependent. However, there is no reason to expect that the higher terms in the cluster expansion 
are negligible. The attractive parts of higher cluster terms can be presented as contributions 
of multi-pion resonances. The properties of the resonances (masses and widths), and even 
their existence themselves, are not very well known. Therefore, an accurate estimation of the 
numerical value of parameter a in the VDW system of pions looks rather problematic, and it 
is treated as purely phenomenological free parameter in the present study. 

V. COMPARISON WITH THE HAGEDORN MODEL 

A concept of the limiting temperature for the hadron gas was introduced in physics by Hage- 
dorn miiiH]. His limiting temperature Tq = 160 ± 10 MeV emerged due to the exponentially 
increasing mass spectrum p{m) for the hadron excited states at large m: 



( 12 ) 


where c, Mq, Tq, and a are the model parameters. These excited states called fireballs were 
considered as a point-like non-interacting particles. A presentation of particle attractive inter¬ 
actions by the resonance states was hrst suggested in Ref. [19] in statistical physics and then 
applied in Ref. 1201 to hadron production. This idea was then further extended within the 
F-matrix formulation of statistical mechanics developed in Ref. 1211 . 

The spectrum p(m) can be found from the statistical bootstrap equation [22] . This equation 
requires the spectrum of low-lying hadron states as an input. In a simplest version of the model, 
this input is reduced to a single lightest hadron - the pion, and the continuous mass spectrum 
p{m) starting from Mq > 2/71,^. The pressure of the Hagedorn model is then defined as 




dmp[m) rn?K 2 {rn/T) , 


(13) 
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where K 2 is the modihed Bessel function. The hrst term in the right hand side of Eq. (13) 
corresponds to the pressure of the ideal Bose gas of pions and the second one to the contribution 


of all excited states (12) taken in the Boltzmann approximation. Other thermodynamical 
functions of the Hagedorn model are then given by the thermodynamical relations ([^. 

From Eqs. (12) and (13) it follows that the thermodynamical functions do not exist at 


T > Tq because of a divergence of the integral with respect to m in Eq. (13). The singular part 
of the integral contributes to the pressure as pn oc (Tq — and to the energy density as 

Eh oc (To — An explicit solution of the statistical bootstrap equation |22] gives a = 3 

[23] . Thus, Eh oc (Tq — and —)■ 00 at T —)■ Tq. On the other hand, for a > 7/2 all 

thermodynamical functions in the Hagedorn model remain finite in a vicinity of the limiting 
temperature T = Tq. Even more, there is a remarkable correspondence between the pion gas 
with the VDW EoS and the Hagedorn model with 7/2 < a < 9/2, i.e., the behavior of pn 
and Eh is very similar to that of p and £ of the VDW EoS shown in the panels (b) and (c) of 
Fig. for a/b = 500 MeV and 600 MeV. In the vicinity of the limiting temperature Tq the heat 
capacity and the speed of sound squared behave in Hagedorn model as (T oc (Tq — and 

cl oc (Tq — respectively. Therefore, for a = 4 the behavior of C and cl in the Hagedorn 

model at T —?■ Tq is the same as in the pion gas with the VDW EoS, if the Tq values in both 
models are set to be equal. In Fig. these two quantities are compared for the VDW model 
with b = 0.45 fm^, a/b = 500 MeV and the Hagedorn model with Tq = 155 MeV and a = 4. 


VI. SUMMARY 

In summary, the interacting pion system has been studied using the quantum version of the 
VDW EoS within the GCE. The role of repulsive and attractive interactions described by the 
VDW parameters b and a, respectively, was considered. It is found that for b = 0.45 fm^ and 
a/b = 200 MeV the VDW EoS for pions appears to be close to that of the ideal pion gas, i.e., 
repulsive and attractive interactions are approximately canceled out. At each value of b there is 
a minimal value of a = a^ such that the limiting temperature T = Tq appears at all a > a^, and 
Tq decreases with increasing a. When T —)■ Tq several remarkable effects happen: the speed of 
sound approaches zero as (Tq — T)^/^, while the heat capacity as well as the scaled variance of 
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Figure 3: The heat capacity (a) and the speed of sound (b). The solid lines correspond the VDW 
EoS with b = 0.45 fm^ and a/b = 500 MeV. The dotted lines correspond to the Hagedorn model ( |l2[ ) 
and (13) with Tq = 155 MeV, a = 4, c = 3.57 fm“^, and Mq = 507 MeV. 


fluctuations of the particle multiplicity go to infinity as (Tq — 

The presented VDW model for pions is compared to the Hagedorn model, where the phe¬ 
nomenon of limiting temperature is also present. Even though the realization of the attractive 
mechanism between particles in the Hagedorn model is different from that in the VDW model, 
we found that a very similar thermodynamical behavior emerges in both models. It takes 
place for the special choice of a power, namely in the pre-exponential factor of the mass 
spectrum p{m) of the Hagedorn model. 

A presence of the limiting temperature Tq in the VDW pion gas is definitely a signal of the 
restricted validity of this model. Similar to the Hagedorn model, the limiting temperature of 
the VDW pion system should be transformed to the temperature of deconflnement transition 
when the fundamental quark-gluon degrees of freedom are introduced. This consideration is, 
however, outside of the scope of the present paper. 
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